Systematics of one-quasiparticle configurations in neutron-rich Sr, Zr, and Mo odd 
isotopes with the Gogny energy density functional. 
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The systematics of one-quasiparticle configurations in neutron-rich Sr, Zr, and Mo odd isotopes is 
studied within the Hartree-Fock-Bogoliubov plus Equal Filling Approximation method preserving 
both axial and time reversal symmetries. Calculations based on the Gogny energy density functional 
with both the standard D1S parametrization and the new DIM incarnation of this functional are 
included in our analysis. The nuclear deformation and shape coexistence inherent to this mass 
region are shown to play a relevant role in the understanding of the spectroscopic features of the 
ground and low-lying one-quasineutron states. 

PACS numbers: 21.60.Jz, 21.10.Pc, 27.60.+j 



I. INTRODUCTION 

The theoretical description of odd-mass nuclei is a dif- 
ficult task in mean field models and, for that reason, they 
have been much less studied than the corresponding even- 
even systems although they roughly constitute half of 
the existing nuclides. The properties of odd nuclei are 
however relevant, for example, to understand odd-even 
effects related to pairing correlations, as well as to con- 
strain nuclear energy density functionals (EDFs) aiming 
to reach a reasonable spectroscopic quality. Within this 
context, the energies, spins and parities of the (neutron 
and/or proton) multi-quasiparticle excitations are basic 
pieces of information to get insight into the underlying 
nuclear structure. Thus, the interest in the properties 
of odd- A nuclei and their computational challenges have 
been lately revived [H-Q. 

One of the main technical difficulties in the descrip- 
tion of odd nuclei is that the (exact) blocking procedure 
0> 0] requires the breaking of time-reversal invariance 
making the calculations troublesome. Moreover, since 
blocking does depend on the choice of the alignment ori- 
entation, for each quasiparticle excitation reorientation 
effects should also be considered d, H| retaining in each 
case the solution with the lowest energy. Therefore, in 
order to reduce the computational cost, in practical appli- 
cations it is convenient to preserve time-reversal as self- 
consistent symmetry. This has been traditionally accom- 
plished within the so called Equal Filling Approximation 
(EFA), according to which the odd nucleon is half into a 
given orbital and half into its time-reversed partner. The 
microscopic justification of the EFA is based on standard 
ideas of quantum statistical mechanics and has been con- 
sidered in Ref. [3J . The predictions arising from various 
treatments of blocking have been compared in Ref. Q. 
On one hand, it has been shown that both exact block- 
ing and EFA approach are strictly equivalent when the 
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time-odd fields of the EDF are neglected. On the other 
hand, it has also been concluded that the EFA is suffi- 
ciently precise for most practical applications. The EFA 
is precisely the approximation employed in the present 
work to study the systematics of the ground and low- 
lying one-quasineutron states in odd- A nuclei with pro- 
ton (neutron) numbers 38 < Z < 42 (47 < N < 67). 

The local systematics of one-quasiparticle excitations 
have already been studied in various mass regions either 
from macroscopic-microscopic approaches [9l4l2l| or from 
the EDF theory d 0101]. To the best of our knowl- 
edge, still the only global study of odd-nuclei within 
the nuclear EDF framework has been reported in Ref. 
[!| where Skyrme-HF+BCS predictions for ground-state 
spins and parities have been compared with finite-range 
droplet model results. 

As already mentioned, in this work we study neutron- 
rich Sr, Zr, and Mo odd isotopes with mass numbers 
A ~ 100. Nuclei in this region of the nuclear chart are 
being actively studied both theoretically [l9l - |22l | and ex- 
perimentally [23M30l ] because of their interest from the 
nuclear structure perspective, as well as because of their 
significance in nuclear astrophysics. In particular, their 
masses and decay properties are an essential input to 
model the path, the isotopic abundances and the time 
scale of the r process in a reliable way [3l| . In addition, 
this region is also characterized by the strong competi- 
tion between various shapes, giving rise to shape insta- 
bilities that lead to shape coexistence and sudden shape 
transitions [2l|, [3^. Thus, it is interes ting to explore the 
predictive power of the Gogny-EDF [33( when applied 
to describe the ground and low-lying one-quasineutron 
states in odd-^4 nuclear systems belonging to this region. 
Although it is clear that one cannot expect a priori to 
reproduce in detail spectroscopic properties of the con- 
sidered odd- A nuclei from global EDFs (see, for example, 
f3~4T - l36j ) the question still remains about the ability of the 
Gogny-type EDFs to describe, at least qualitatively, the 
main features of the observed trends in the Sr, Zr and 
Mo isotopic chains. 
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In addition to D1S which is the most standard 
and thoroughly tested parameter set (see, for example, 
[38T - I53T ] and references therein), in the present study we 
have also considered the most recent parametrization of 
the Gogny-EDF, namely the DIM [54[ parametrization. 
The goal is twofold, first we want to verify the robustness 
of our predictions with respect to the particular version 
of the EDF employed. If the results with two different 
EDFs (or two parametrizations of one, as is our case) 
are similar, we can be rather confident about the inde- 
pendence of the predictions with respect to the details 
of the EDF. The second goal is to test the performance 
of DIM in the present context (spectroscopy of odd-A 
nuclei) as to add arguments to decide whether DIM [HiJ 
can supersede the rather old D1S (37j parametrization 
or not. The fitting protocol of DIM includes input from 
a realistic Equation of State (EoS) both in symmetric 
and neutron nuclear matter. It also explores, in the sym- 
metric nuclear matter EoS case, the four possible spin- 
isospin channels and tries to reproduce the trends ob- 
tained in realistic Brueckner-Hartree-Fock calculations. 
As a consequence, the DIM parametrization reproduces 
the correct sign for the isovector splitting of the effective 
mass for neutron-rich matter at all possible asymmetries. 
The consequence of such a sophisticated fitting protocol 
is an energy root mean square deviation for the exper- 
imentally known 2149 masses of only 0.8 MeV which is 
much better than the one of D1S. The systematic com- 
parison [HH of the collective 2 + excitation energies ob- 
tained with D1S and DIM (only even-even nuclei were 
considered) suggests that the later essentially keeps the 
same nuclear structure predictive power of the former. 
Subsequent analysis centered on the structure of the in- 
trinsic wave functions [55| has confirmed the previous 
findings. Therefore, it is timely and necessary to fur- 
ther extend the comparison between D1S and DIM to 
odd nuclei. To this end, we analyze in the present pa- 
per the spectroscopic properties of neutron-rich Sr, Zr 
and Mo odd isotopes in the A ~ 100 mass region within 
the Hartree-Fock-Bogoliubov (HFB) Q plus EFA frame- 
work [3| (HFB-EFA). The election of the considered nu- 
clei is mainly driven by the intense experimental (23l - [29l | 
and theoretical |19l422l | efforts to better characterize the 
structural evolution of the ground and excited nuclear 
shapes in this region of the nuclear chart [56| . 

The paper is organized as follows. In Sec. HH we 
present a brief description of the theoretical formalism 
used in the present work, i.e., the HFB-EFA framework. 
The results of our calculations for the considered nuclei 
are discussed in Sec. Mil where we pay attention to mean 
field potential energy curves (PECs), equilibrium defor- 
mations of the various competing shapes, single-particle 
energies, one-quasineutron states and their spectroscopic 
evolution along the Sr, Zr and Mo isotopic chains. We 
also compare our results with the available experimen- 
tal data. Finally, Sec. IIVI is devoted to the concluding 
remarks and work perspectives. 



II. THEORETICAL FRAMEWORK 

In the present study several Sr, Zr and Mo isotopes 
with odd neutron number are studied. Nevertheless, the 
blocking procedure changing the number parity @, @] of 
a given even-even HFB vacuum and providing us with 
the corresponding ground and low-lying one-quasiparticle 
states in odd-A nuclei also requires to consider even-A 
nuclei in these isotopic chains, as discussed below. 

In previous studies of even-even nuclei we have found 
advanta geou s to use the so called gradient method 
[55l [57l |58| to obtain the solution of the HFB equa- 
tions, leading to the (even number parity) vacuum |$). 
Within this method, the HFB equation is recast in terms 
of a minimization (variational) process of the mean field 
energy and the Thouless parameters defining the most 
general HFB wave functions Q are used as variational 
parameters. For a given point in the multidimensional 
variational space of the Thouless parameters the proce- 
dure uses the direction of the gradient in its search for the 
minimum. The advantage of the gradient method over 
others, like the successive iteration method, relies on the 
way constraints are implemented, which allows a larger 
number of them to be treated at once. As it is customary 
in calculations with the Gogny force, the kinetic energy 
of the center of mass motion has been subtracted from 
the Routhian to be minimized in order to ensure that the 
center of mass is kept at rest. The exchange Coulomb en- 
ergy was considered in the Slater approximation and we 
neglected the contribution of the Coulomb interaction to 
the pairing field. Both axial and time-reversal are self- 
consistent symmetries in our calculations. 

The quasiparticle creation and annihilation operators 
k and Pk) associated with a given (even-even) HFB 
vacuum |$) have been expanded in an axially symmetric 
harmonic oscillator (HO) basis (c,,C() 

tn 

= y ] U mfl c m + V m ^c^ n . (1) 

tn 

containing enough number of shells (13 major shells) to 
grant convergence for all values considered of the mass 
quadrupole moment Q20 = ^(^\2z 2 -x 2 — y 2 \<&) [58[ and 
for all the nuclei studied. 

Constrained calculations have been performed to gen- 
erate PECs for even-even Sr, Zr and Mo nuclei. Such 
PECs arc displayed in Fig. Q] for a set of representa- 
tive isotopes and the parametrization Gogny-DlS. One 
can also find a systematic compilation of PECs obtained 
with Gogny-DlS in Ref. (5j|. The computation of such 
PECs is twofold: first, they give us initial hints on the 
evolution of the different competing shapes in the con- 
sidered nuclei and, second, they provide a whole set of 
prolate, spherical and oblate (reference) even-even HFB 
states for the subsequent treatment of the neighboring 
odd- TV nuclei. 
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FIG. 1. (Color online) PECs for Sr (a), Zr (b), and Mo (c) 
isotopes with N =52, 56, 60, 64, and 68 obtained from HFB 
calculations with the Gogny-DIS energy density functional. 



In fact, once a reference (even-even) state with a given 
deformation is chosen, we use it to perform an additional 
(constrained) HFB calculation providing an unblocked 
fully-paired state corresponding to an odd average neu- 
tron number (false vacuum with the same defor- 
mation. Such prolate, spherical and oblate false vacua 
are then used as input configurations in our subsequent 
blocking scheme (i.e., EFA). 

The HFB ground-state wave function |4>) of an even- 
even nucleus is defined by the condition of being the vac- 
uum of the annihilation quasiparticle operators (3^ de- 
fined in Eq. (|T|). Regarding observables (mean values) 
the HFB ground state of an even-even system is specified 
7] by the (hermitian) density matrix 



and the (antisymmetric) pairing tensor 
KU = (V*U T ) :1 



(2) 



(3) 



where the U and V amplitudes are those of Eq ([T]) defin- 
ing the Bogoliubov transformation. Applying any uni- 
tary transformation to the right of both amplitudes U 
and V leaves both p and n (and therefore any mean 
value) unaltered whereas the quasiparticle operators are 
affected by the representation of the unitary transforma- 
tion in operator space. This invariance is not maintained 
in the odd- A case, as we will discuss below, implying con- 
ceptual changes in the treatment of an odd- A system as 
compared to the simpler even-even case. 

On the other hand, the ground and low-lying one- 
quasiparticle states of odd- A systems, like the ones con- 
sidered in the present work, can be handled with blocked 
(odd number parity [1,0]) HFB wave functions 



(4) 



where fis indicates the quasiparticle state to be blocked 
and stands for the indexes compatible with the symme- 
tries of the odd-nuclei, such as the angular momentum 
projection K and parity in the case of axial symmetry. 
In practice, the blocked state (Q| can be accomplished 
by exchanging the columns labeled by the index ps in 
the HFB amplitudes U and V [see, Eq. (JTJ] with the 
same columns in the matrices V* and U*, respectively 

The density matrix and the pairing tensor correspond- 
ing to the one-quasiparticle state @, read 



and 



.(Mi 



V*V 2 



V^ 1 



TT- V* 



(5) 



(6) 



They violate time-reversal invariance and so do the HF 
and pairing fields associated with them, making the nu- 
merical calculation much more difficult than in the even- 
even case. These matrices are no longer invariant under 
unitary transformations of the U and V amplitudes ap- 
plied to their right. As a consequence, that degree of 
freedom has to be taken into account in the solution of 
the HFB equations. A manifestation of this effect is the 
reorientation effect [4|, |8[ discussed recently. 

A very useful approximation to (exact) blocking pre- 
serving time-reversal invariance is the EFA Hi], where 
one writes the density matrix 



(EFA,hb) 
Hij 



V* V- 



+ l( U ^ B U h B ~ V ^ B ) (7) 
and the pairing tensor 
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FIG. 2. (Color online) Single particle energies for protons 
and neutrons in 102 Zr as a function of the axial quadrupole 
moment Q20- The Fermi level is depicted as a thick dashed 
red line. The results have been obtained with the Gogny- 
D1S EDF. Solid lines correspond to levels with positive 
parity whereas dashed lines correspond to negative parity 
states. Asymptotic (Nilsson) quantum numbers [N,n z ,A]K w 
are shown in the neutron case for Q20 values close to those 
where the minima of the PECs are located (vertical arrows). 



l(um B v^ B -v± B U J1Zl 



(8) 



including not only the [Ib state but also its time reversed 
partner ~p, B . The corresponding HF and pairing fields 



/ j V i<13pPpq 
qp 



(EFA^ b ) _ 1 \p (EFA^ B ) 
qp 



(9) 



(10) 



and the total EFA energy can be written in the usual 
HFB form 



E (EFA^ B ) _ Tr 

+ \Tr 

-'-Tr 
2 



tp (EFA,p B 

Y(EFA,pb) p (EFA,p, B ) 
^(EFA,fi B ) k (EFA,hb) 



(11) 



The microscopic justification of this expression for 

E (EFA,vb) has been first 

given in Ref. Q using ideas 
of quantum statistical mechanics. The EFA energy can 
be obtained as the statistical average, with a given den- 
sity matrix operator, and applying the variational prin- 
ciple to it, the usual HFB-EFA equation Q is obtained. 
The existence of a variational principle allows to use the 
gradient method to solve the HFB-EFA equation [3[ with 
the subsequent simplification in the treatment of the con- 
straints. 

The solution of the HFB-EFA (as well as the exact 
blocked HFB) equation depends upon the initial blocked 
level [Ib ■ In the HFB-EFA case, the K quantum number 
is selfconsistently preserved along the calculation and so 
does the parity if octupole correlations are not allowed 
in the iterative process. Blocking levels with different K 
values and parities lead to different quantum states of 
the odd- A nucleus being the ground state the one with 
the lowest energy. Given a K value and parity, one could 
naively think that considering the quasiparticle with the 
lowest quasiparticle energy as the initial blocked level 
should lead to the lowest energy solution for that value of 
K and parity. However, due to the selfconsistent nature 
of the whole process this is not by any means guaranteed 
and therefore several quasiparticles, usually correspond- 
ing to the lowest quasiparticle energies have to be con- 
sidered. In addition, in the present case and owing to 
the presence of coexisting prolate, oblate and spherical 
minima in some of the nuclei considered, blocked config- 
urations with those quadrupole deformations have to be 
explored. For this reason the minimization process has 
to be carried out several times using different initial pro- 
late, spherical and oblate (false) vacua. We have repeated 
each calculation, for a given false vacuum and JsT-values 
from 1/2 up to 15/2, using as initial blocking states the 
12 quasiparticles corresponding to the lowest quasiparti- 
cle energies. The use of so many initial configurations is 
to make sure we are not going to miss the true ground 
state and all the lowest excited states. The resulting com- 
putational cost is high as in the worst case (three false 
vacua) we have to carry out for each nucleus a total of 
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FIG. 3. (Color online) Experimental (a) excitation energies 
and spin-parity assignments of the non-collective states are 
compared with HFB-EFA results (b) for the one-quasineutron 
states in odd-iV Sr isotopes (see text for details). Prolate 
configurations are shown by black lines, oblate ones by red 
lines, and spherical ones by blue lines. 



8 x 3 x 12 = 288 HFB-EFA calculations. Among all the re- 
sulting one-quasineutron configurations, we have selected 
those within a 1 MeV energy window and compare them 
with available experimental data (see, for example, Figs. 
El, SI and [S] below). Note that for nuclei in this region of 
the nuclear chart, there exist several competing shapes 
at low excitation energy and therefore our procedure as- 
sures that the lowest energy solution can be reached for 
all values of the quadrupole moment Q20 an d mass num- 
ber. Finally, let us also mention that in those cases where 
the neutron pairing energy of a given false vacuum is too 
small, we increase it by constraining the mean value of 
the number particle fluctuation neutron operator (AN) 
in this false vacuum. We have found that this procedure 
leads to a very fast convergence of the calculations. 



III. RESULTS 

A. Potential energy curves and single-particle 
energies 

It has been shown experimentally that the ground 
states of even-even Sr, Zr and Mo isotopes with N rang- 
ing from the magic neutron number N = 50 up to N ~ 60 
are weakly deformed, but they undergo a shape transition 
from nearly spherical to well deformed prolate (oblate) 
configurations as N = 60 is approached and crossed. Evi- 
dence for such an abrupt shape /phase transition includes 
2 + lifetime measurements 1241. |28[ and quadrupole mo- 
ments for rotational bands 1251] , as well as isotopic shifts 
in nuclear charge radii (23l. l26l |29| . Heavier Sr and Zr 
(A ~ 110) isotopes exhibit axially symmetric well de- 
formed shapes, whereas the heavier Mo isotopes display 
a tendency towards triaxiality [2l],|27|- Above this region, 
it has been suggested 6l| that the N = 82 shell closure 
might be quenched far from stability but still weak ex- 
perimental evidence has been found. 

In Fig. [TJ we show the PECs for the even-even Sr 
(a), Zr (b), and Mo (c) isotopes. We display PECs cor- 
responding to N =52, 56, 60, 64 and 68 isotones from 
which, the main features of the shape evolution can be 
followed. Nuclei with N = 52 show a sharp spherical 
minimum that becomes rather shallow at N — 56. Iso- 
topes with N — 60 are already deformed with oblate 
and prolate minima very close in energy. In the case of 
Sr isotopes the ground state is prolate, for Zr isotopes 
oblate and prolate minima are almost degenerate, while 
the ground state becomes oblate for Mo nuclei. For neu- 
tron numbers greater than 60, prolate and oblate min- 
ima become well defined. Our results agree well with 
the shape/phase change predicted around N = 60 by 
microscopic-macroscopic models fl 9| . [2 fj|. I62L [63j and mi- 
croscopic selfconsistent relativistic [64[ and non relativis- 
tic [13 IsH mean field calculations. Later on, in section 
IIIIB1 our analysis of the shape evolution in even-even Sr, 
Zr and Mo nuclei will be further extended with the sys- 
tematics of the ground and low-lying one-quasineutron 
configurations predicted by our HFB-EFA calculations 
for odd- A nuclei in these isotopic chains. 

Another interesting outcome of our HFB calculations 
for even-even nuclei (see, also the discussion in section 
IIII B[) . are the proton and neutron single-particle ener- 
gies (SPEs) shown in Fig. [2] as functions of the axial 
quadrupole moment Q20, for the case of 102 Zr (Z =40, 
N =62). To obtain these Nilsson-like diagrams, we 
have computed the eigenvalues of the Routhian h = 
t + r — A20Q20, with t being the kinetic energy operator, 
r the Hartree-Fock field and A20Q20 the term containing 
the Lagrange multiplier used to enforce the correspond- 
ing quadrupole constraint. Obviously, the usual mean 
field constraints Q on both neutron and proton numbers 
are also taken into account. Fermi levels are plotted with 
thick dashed (red) lines in the figure. These type of di- 
agrams help us [55|, HH to identify regions of low level 
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FIG. 4. (Color online) Same as in Fig. but for odd-TV Zr 
isotopes. 



density around the corresponding Fermi levels which, ac- 
cording to the Jahn- Teller effect [66[, favor the onset of 
deformation. 

For axially symmetric configurations, the SPE levels 
are tagged by the (half integer) K quantum number that 
corresponds to the third component of the angular mo- 
mentum in the intrinsic frame. Due to the time-reversal 
invariance, single-particle orbitals with the same absolute 
K value are degenerate (Kramers degeneracy). Positive 
(negative) parity levels are plotted with solid (dashed) 
lines. At the spherical configuration Q20 = 0, the quan- 
tum numbers nlj are recovered. The asymptotic quan- 
tum numbers [N, n z , A]K n are shown in Fig. [2] for Q20 
values close to those where the minima of the PECs (see, 
Fig. [TJ are located (vertical arrows). For each single 
particle state they have been assigned according to the 
largest amplitude in the expansion of the single particle 
wave function in the axially symmetric harmonic oscilla- 
tor basis. 

For the neutron-rich A ~ 100 nuclei of our concern in 
the present work, the valence protons occupy the N = 3 
shell and start to fill the g 9 / 2 orbitals in the N = 4 
shell. Neutrons occupy states belonging to the N = 4 
shell, which are strongly mixed by deformation, and be- 
gin to fill the ft.11/2 intruder orbitals coming down from 
the N = 5 oscillator shell. Because of their opposite par- 
ity, the h±i/2 neutron states (the gg/ 2 proton states) do 



not mix with the N — 4 neutron orbitals (the N = 3 
proton orbitals). These intruder orbitals polarize the 
core towards oblate and prolate deformations. As a con- 
sequence, the underlying nuclear structure in this mass 
region is very sensitive to the occupancy of these single- 
particle orbitals and the result is a rapid change in nu- 
clear spectroscopic properties as a function of both neu- 
tron and proton numbers. 

Below we compare the results obtained with the D1S 
parametrization of the force with those obtained with the 
newly introduced DIM parametrization. As it turns out 
that the spin orbit strength of DIM is 12 % smaller than 
the one of D1S one may wonder whether this difference 
could lead to a strong impact in the single particle spec- 
trum and its behavior with deformation. To answer this 
question we have plotted the DIM single particle levels 
and compared them to the D1S ones. The comparison 
of the two sets of levels shows small differences in the 
position of the levels (of the order of a few hundred of 
keV at most) that are not strong enough as to change the 
ordering of the spherical single particle spectrum. The 
behavior with deformation is not changed in a significant 
way either. This small impact can be related to the larger 
effective mass of DIM as compared to D1S, that makes 
the single particle spectrum more dense in DIM than in 
D1S. As a consequence a weaker spin-orbit strength is 
required in DIM to shift the intruders to the lower shell. 



One-quasiparticle states in odd-iV Sr, Zr and 
Mo nuclei 



In Figs. [3j |4] and [5] results are given for odd-A^ Sr, 
Zr, and Mo isotopes, respectively. In each figure the ex- 
perimental [56| ground and low-lying one-quasineutron 
states, characterized by their spin-parity assignments 
(panels (a) of the figures) are compared with the ones 
predicted within our HFB-EFA calculations (panels (b) 
of the figures) . This is certainly a very challenging test 
for the spectroscopic quality of the Gogny-DlS EDF for 
which results are shown. It is worth noting that in 
some cases the experimental spin and parity assignments 
in panels (a) have been deduced from systematics and 
should be treated with caution. 

The existence of oblate and prolate competing shapes 
in this mass region prevents us to use the usual conven- 
tion [63| of plotting the theoretical results with hole 
(particle) states below (above) the corresponding ground 
state. In Figs. [31 S and [SJ the excited states in a given 
isotope are referred to the corresponding ground state re- 
gardless its shape. Prolate configurations are shown by 
black lines, oblate ones by red lines, and spherical ones 
by blue lines. The quasiparticle states in the deformed 
configurations are labeled by their K n quantum numbers 
and, in addition, the ground states are labeled by their 
asymptotic (Nilsson) quantum numbers [N, n Zl A] AT 71 ". In 
the lighter isotopes we obtain spherical (sharp or shal- 
low) equilibrium shapes and therefore only the spherical 
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orbital, which is either gg/ 2 or d$/2 in all the cases, is in- 
dicated. Like-orbitals in neighbor isotopes are connected 
by dashed lines to better appreciate their isotopic evolu- 
tion. 

Looking at Figs. |4] and the striking correlation 
between the ground state deformations in the even-even 
isotopes (Fig. [IJ and the ground state shapes of the odd 
ones is observed. The lighter Sr nuclei (A = 85 — 91) are 
spherical. They become oblate for A = 93 and 95 with 
prolate states at low excitation energies. Finally, prolate 
deformed ground states are obtained in the heavier iso- 
topes with A — 97 — 105. Oblate excited states below 1 
MeV are obtained for A = 97, whereas they lie above 1 
MeV in the heavier ones. 

The lighter Zr isotopes (A = 87 — 95) displayed in 
Fig. |4l show spherical equilibrium shapes, which turn 
into oblate (A = 97,99), prolate (A = 101 - 105), and 
again oblate (A = 107) ground states with excited one- 
quasiparticle configurations of different shapes below 1 
MeV and very close in energy, in agreement with Fig. Q] 
(b). 

The odd Mo isotopes in Fig. [5] also display similarities 
with Fig. [1] (c) . Again, spherical ground states are found 
for the lighter isotopes. Alternating oblate and prolate 
equilibrium shapes are predicted in 99 ! 101 > 103 > 105 Mo. Pro- 
late and oblate configurations are particularly close in 
energy within this mass range. Oblate ground state de- 
formations are found for both 107 > 109 Mo with no prolate 
one-quasineutron configurations in 109 Mo below 1 MeV. 

In Fig. [3l we compare the experimental (a) low-lying 
one-quasineutron spectra with the ones predicted by our 
HFB-EFA calculations (b) for odd- TV Sr isotopes. We can 
see a reasonable agreement in the ground states, except 
for the isotopes 93 > 95 > 97 Sr. In 93 Sr, the 5/2 + experimen- 
tal ground state appears at 0.1 MeV in our calculations, 
whereas in 95 - 97 Sr the l/2 + experimental ground states 
are found as oblate excited configurations at E = 0.25 
and E = 1 MeV, respectively. The ground states of the 
heavier isotopes "• 101 Sr are well reproduced. 

The structural evolution of the ground and excited one- 
quasiparticle states along the considered Sr, Zr and Mo 
isotopic chains can also be followed by looking at the 
dashed lines connecting the states characterized by the 
same asymptotic quantum numbers [N, n z , A]K n and at 
Fig. [21 as a guide. Indeed, the ground state sequence as 
the number of neutron increases exhibits a clear corre- 
spondence with the orbitals found in Fig. [2] along thick 
vertical arrows at oblate (Q20 ~ —2.5 b) and prolate 
(Q20 ~ 5 b) quadrupole moments. For example, one 
sees how the ground state in 101 Sr, which corresponds 
to the prolate [532] 5/2" state, becomes the first one- 
quasineutron excitation, below 0.6 MeV, in the neighbor- 
ing 99 Sr and 103 Sr with prolate [411]3/2+ and [411]l/2+ 
ground states. On the other hand, 97 Sr exhibits a prolate 
[404] 9/ 2 + ground state, which turns out to be a quasineu- 
tron band head in 95 Sr and 99 Sr while the [411]l/2 + pro- 
late ground state in 103 Sr appears as an excited state 
between 0.6 and 0.4 MeV in 101 Sr and 105 Sr, respectively 



In Fig. [4] the comparison is made for Zr isotopes. The 
agreement is again very reasonable exception made of 
97,99 Zr for which the reported ground states appear in 
our HFB-EFA calculations as oblate excited states be- 
low 0.5 MeV. Special mention deserves, the remarkable 
agreement for the heavier isotopes 101 ' 103 Zr, where our 
calculations reproduce both the ground states and the 
low-lying excitations. One can see, how the band head 
5/2" at 0.45 MeV in 101 Zr and the ground state in 103 Zr 
correspond to the same [532]5/2~ single-particle config- 
uration, as it has been experimentally established in Ref. 
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FIG. 5. (Color online) Same as in Fig. O but for odd-TV Mo 
isotopes. 

The dashed lines in FigJH are also useful to trace the 
isotopic evolution of the single-particle configurations, 
which can be also understood by moving the (neutron) 
Fermi level in Fig. [2] through the vertical arrows. Thus, 
taking for example the prolate solutions, if one looks at 
the orbitals intersecting the vertical arrow at Q20 = 5 b 
in Fig. [2l one finds the deepest state at about e sp = —9 
MeV, which corresponds to a [420]l/2 + (si/ 2 ) state. It 
can be seen, in FigJH as the prolate excited state at 
0.6 MeV in 95 Zr. The next prolate states in Fig. [H 
are [550]l/2 _ falling down very quickly from hu/2 and 
[422] 3/2+ (^5/2) appearing as the first prolate excitations 
in 97 Zr. On the other hand, in 99 Zr the first two prolate 
states 3/2~ and 9/2 + appear around 0.2 MeV. They are 
associated with the [541]3/2 _ coming down from h\\/2 
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and [404]9/2+ raising from g g / 2 in Fig. [2] Going fur- 
ther up in Fig. |2 we find [411] 3/2+ (g 7 / 2 ), which is the 
ground state in 101 Zr and [532]5/2 _ (^11/2) which rep- 
resents an excited state at about 0.45 MeV. Further up 
in excitation energy one finds, the 9/2+ and 3/2~ states 
already discussed (see the connections with dashed lines 
with their partners in 99 Zr). In 103 Zr the 3/2+ and 5/2~ 
states interchange their positions with respect to 101 Zr 
and, as it can be easily understood from Fig. [21 the 
[532]5/2~ configuration becomes now the ground state. 
Subsequently, the [411] 1/2+ (57/2) ground state is pre- 
dicted by our HFB-EFA. The excited one-quasineutron 
5/2~ and 5/2 + states, at about 0.35 MeV, correspond to 
the already discussed [532]5/2~ (/in/2), which now lies 
below the Fermi level, and to the [413] 5/2+ (d 5 / 2 ) config- 
uration, above the Fermi level, that becomes the lowest 
prolate deformed one-quasiparticle excitation in the last 
odd- TV nucleus considered (i.e., 107 Zr). 

A similar analysis can be made for the oblate one- 
quasineutron states by looking at the intersections of the 
single-particle levels in Fig. [2] with the vertical line at 
Q20 = —2.5 b. As can be seen from Fig. @J 97 Zr dis- 
plays a ground 3/2+ and two low- lying excited (1/2+ 
and 7/2+) oblate deformed states. They can be associ- 
ated with the three orbitals at -8 MeV in the oblate sector 
of Fig. m namely, [411]3/2+ (d 5/2 ), [411]l/2+ (d B / 2 ), and 
[413J7/2+ {g 7/2 ). Next, in "Zr we find an oblate 3/2+ 
ground state and the oblate 11/2" excitation at 0.1 MeV 
which correspond to the [422]3/2+ (g 7/2 ) and [505] 11/2" 
intruder state coming down from the hn/ 2 shell. The 
[420] 1/2+ and [422]5/2+ (coming from g 7/2 shell) exci- 
tations , are also visible in the quasiparticle spectrum of 
99 Zr. The nucleus 101 Zr exhibits a prolate ground state 
but one finds a [505]ll/2~ as the lowest oblate quasineu- 
tron configuration very close to the 3/2 + excitation (that 
was the ground state in 99 Zr). The same comments ap- 
ply to both 103 Zr and 105 Zr but now the state [420] 1/2+ 
(si/2) comes into play. Finally, 107 Zr displays an oblate 
ground state that can be associated with a [505]9/2~ 
configuration coming down from the ft.11/2 shell. 

According to the discussions in Ref. 



99 c 



101, 



2JJ, where the 



rotational bands of the lowest quasiparticle states were 
measured, and the results of our recent study in Ref. [2l[ , 
where triaxial shapes were predicted in neutron-rich Mo 
isotopes, the odd Mo nuclei are the most difficult to ana- 
lyze. Our HFB-EFA description in terms of axial shapes 
nicely reproduces the ground 5/2 + and the two excited 
3/2+ and 1/2+ states in 97 Mo. Our calculations also pro- 
vide us with the observed 1/2+, 5/2+, 3/2+, and 11/2" 
one-quasiparticle configurations in 99 Mo although redis- 
tributed. The observed 1/2+ , 5/2+, and 11/2" states 
are obtained in 101 Mo within 0.4 MeV excitation energy. 
On the other hand, the ground state in 103 Mo is well 
reproduced while the excited 5/2 _ appears at about 0.6 
MeV. The latter becomes the experimental ground state 
in 105 Mo but it is predicted at 0.15 MeV in our calcula- 
tions. Finally, the observed ground states in both 107 Mo 
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FIG. 6. (Color online) Excitation spectra of some Sr isotopes 
obtained within the HFB-EFA calculations based on the D1S 
and DIM parametrizations of the Gogny-EDF are compared 
with available experimental data. The color scheme used to 
characterize the quadrupole deformation of each of the states 
is the same as in Fig. [3] Experimental levels are plotted in 
green because there is no clear experimental evidence con- 
cerning the sign of their intrinsic deformations. 



and 109 Mo, are obtained as oblate excited (about 0.5 
MeV) configurations within our HFB-EFA calculations. 

In general, we observe that the main difficulties to 
describe the experimental information appear in those 
cases where the mean field approach might not be suffi- 
cient. Such is the case of N — 55 — 59 isotopes, where 
we get very shallow minima (see Fig. [I}, a situation 
that in general requires configuration mixing for a better 
description. Work along this direction has already been 
carried out in [68, 69] in the context of a collective hamil- 
tonian obtained in the Gaussian Overlap Approximation 
(GOA) framework. We believe that a better treatment 
of the problem is required to pin down the subtle details 
of the spectrum of an odd- A system and therefore the ex- 
act Generator Coordinate Method (GCM) with the axial 
quadrupole moment would be required. Work along this 
line is in progress and will be reported elsewhere. Dif- 
ficulties also appear when the description requires the 
consideration of triaxility as it is the case for the heavier 
Mo isotopes, which exhibit a 7-soft behavior and even 
triaxial minima [2l| with the axial ones converted into 
saddle points. In such a situation one also expects that 
the oblate and prolate solutions will be highly mixed with 
the intermediate triaxial configurations. On the other 
hand, a HFB axial description is expected to work bet- 
ter for well developed axial minima separated by energy 
barriers in both spherical and triaxial shapes as in the 
heavier isotopes of Sr and Zr. 

Another interesting point worth noticing is the isotonic 
behavior of the considered odd- A nuclei. Because the iso- 
topes studied in this work have an odd number of neu- 
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FIG. 7. (Color online) Same as in Fig. [6] but for Zr isotopes. 



trons, one observes very similar spectroscopic schemes 
for fixed number of neutrons (isotones), where the odd 
neutron determines the spin and parity of the ground 
state. First, one realizes the observed experimental cor- 
respondence between the spins and parities of isotones 
within Sr, Zr, and Mo nuclei. Thus, starting at N = 47 
with ,r = 9/2+ in 85 Sr, 87 Zr, and 89 Mo, the sequence 
of 9/2+, 9/2+, 5/2+, 5/2+, 5/2+, 1/2+, 1/2+, 3/2+, 5/2" 
is observed experimentally from TV = 47 up to N = 63. 
The agreement between these experimental findings and 
our calculations is specially remarkable in the heavier iso- 
tones. 

For N = 61 ("Sr, 101 Zr, 103 Mo) the ground state cor- 
responds to [411] 3/2+ states both experimentally and 
theoretically. In addition, a low- lying band head 5/2~ 
is observed experimentally |27J in 101 Zr and 103 Mo that 
corresponds reasonably well with the predictions of our 
HFB-EFA calculations. The experimental assignment 
(5/2+) to the excited state in 99 Sr is uncertain and could 
also correspond to a 5/2~ in which case the similarity 
would be complete. The N = 63 isotones 101 Sr and 
103 Zr exhibit a [532]5/2~ ground state in agreement with 
experiment. In the case of 105 Mo we obtain an oblate 
ground state [420]l/2+, which does not show up ex- 
perimentally, but the first excited prolate state is again 
[532] 5/2 ~. Excited 3/2+ quasiparticle configurations ob- 
served experimentally at N = 63 are also well repro- 
duced in the calculations as prolate configurations. The 
N = 65 isotones 103 Sr and 105 Zr are predicted to be 
prolate [411] 1/2+ in their ground states, while 107 Mo is 
predicted as an oblate [420] 1/2+ at variance with exper- 
iment. 

The main difficulty to understand the previous corre- 
spondence appears in the isotones N =57, 59, where the 
experimental 1/2+ assignments to their ground states are 
difficult to reproduce in our calculations. However, these 
1/2+ quasiparticle states appear in the calculations as 
excited oblate configurations at low excitation energies. 



The 3/2~ excited state observed experimentally in 97 Sr 
and 99 Zr appear in our calculations at similar energies 
as prolate configurations. In 99 ~ 101 Mo there are two ex- 
cited 5/2+ and 11/2" states also found in our HFB-EFA 
calculations below the 1 MeV energy window of the plot. 

As mentioned in Sec. HI one of our aims in the present 
study was to further extend the comparison, already un- 
dertaken in Ref. [55|, between the performance of the 
parametrizations Gogny-DlS and Gogny-DIM to the de- 
scription of the spectroscopy of odd-^4 nuclei. There- 
fore, in Fig. ©(for 99 ' 101 Sr), Fig. [7] (for 101 ' 103 Zr) and 
Fig. [H] (for 103 ' 105 Mo), we compare the results obtained 
for the low-lying one-quasineutron states using these two 
different incarnations of the Gogny-EDF with the exper- 
imental data. First, one realizes the nice similarity be- 
tween the HFB-EFA predictions obtained with these two 
parametrizations. In particular, the predicted ground 
states coincide in most of the cases and agree well with 
experiment. Only in Mo isotopes one observes some dis- 
agreement but in this case, we obtain a high density of 
very low-lying configurations with the observed ground 
states always below 0.2 MeV. 
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FIG. 8. (Color online) Same as in Fig. [6] but for Mo isotopes. 

Looking at the figures in more detail, we can see that 
within the considered energy window of 0.7 MeV, practi- 
cally the same states are obtained from the calculations 
with D1S and DIM. This is especially true for prolate 
deformed states. For example, in Sr isotopes (Fig. [6]) we 
obtain the same 5/2 ~ states in 99 Sr as well as the 3/2+ 
and 1/2+ states in 101 Sr very close to the experiment and 
with very small variations (less than 0.2 MeV) from one 
EDF to another. 

The same happens in Fig. [7] for the Zr isotopes. The 
prolate states 5/2" (in 101 Zr) and 3/2+, 1/2+ (in 103 Zr) 
are predicted very close with D1S and DIM and agree 
well with experiment. However, we observe an enhanced 
sensitivity in the case of the oblate configurations, result- 
ing in lower energies, of the order of a few hundred keV, 
in the case of DIM with respect to D1S. This is a conse- 
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quence of having closer oblate and prolate minima with 
the parametrization DIM. The ground states are nicely 
reproduced. 

In the case of Mo isotopes (Fig. [5J the experimental 
ground and first excited state are interchanged in 103 Mo 
and 105 Mo. Our HFB-EFA calculations provide prolate 
3/2+ and 5/2~ states compatible with experiment, but 
also other oblate configurations at very low excitation 
energy. We can see 11/2", 3/2+, 5/2+, and 1/2+ low- 
lying oblate states in 103 Mo. In addition to those we 
also find an oblate 9/2~ state in 105 Mo. As for Zr iso- 
topes, we also find that DIM tends to favor oblate so- 
lutions, which appear displaced to lower energies with 
respect to the prolate solutions in D IS. In general we ob- 
tain a comparable spectroscopic quality of the two Gogny 
parametrizations and it is very satisfying to observe how 
the new parametrization DIM, in spite of the relaxation 
of some of the original constraints in its fitting protocol 
and more oriented to reproducing nuclear masses [5~3 |. 
still follows very closely the fine details predicted with 
Gogny-DlS for odd-yl nuclei in a region of the nuclear 
chart with such a challenging shape evolution. 

IV. CONCLUSIONS 

In this work we have studied the systematics of one- 
quasineutron configurations in odd-A Sr, Zr and Mo iso- 
topes within the selfconsistent HFB plus EFA frame- 
work, an approach that has been shown to be fairly ad- 
equate for most purposes. However, we are aware of the 
challenge of reproducing in detail the observed spectro- 
scopic properties in the particular mass region consid- 
ered in the present study. Therefore, our aim has been 
to understand qualitatively the structural evolution of 
the ground and low- lying one quasiparticle configurations 



with neutron number. We have also shown that the qual- 
ity of the spectroscopic results obtained with the recent 
DIM parametrization of the Gogny force is comparable 
to the one obtained with the standard parametrization 
D1S. We conclude that both DIM and D1S parametriza- 
tions reproduce, at least qualitatively, the main features 
observed in the isotopic (and isotonic) trends of the 
neutron-rich and odd- A Sr, Zr and Mo nuclei. 

The main deficiency of our mean field description is 
the preservation of axial symmetry. We have found this 
to be relevant in nuclei characterized by shallow minima 
around the spherical configurations like in the A^=55-59 
isotopes, as well as in those nuclei with some tendency to 
7-soft behavior or even triaxial minima like in the heavier 
Mo isotopes [2l[. A beyond mean- field treatment seems 
to be necessary to improve the quality of the description 
in those cases, and configuration mixing calculations in 
the spirit of the Generator-Coordinate-Method (GCM) 
method |7| may be required to improve the mean-field 
results. We are still far from being able to apply such a 
configuration mixing approach within an exact blocking 
scheme for odd nuclei. Nevertheless, a HFB-EFA-GCM 
scheme can be considered as a plausible step forward in 
this direction and work along this lines is in progress. 
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